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Abstract 

We present an accurate investigation of the algebraic conditions that the 
symbols of a convergent, univariate, binary, non-stationary subdivision scheme 
should fulfill in order to reproduce spaces of exponential polynomials. A sub- 
division scheme is said to possess the property of reproducing exponential 
polynomials if, for any initial data uniformly sampled from some exponential 
polynomial function, the scheme yields the same function in the limit. The 
importance of this property is due to the fact that several functions obtained 
as combinations of exponential polynomials (such as conic sections, spirals or 
special trigonometric and hyperbolic functions) are of great interest in graph- 
ical and engineering applications. Since the space of exponential polynomials 
trivially includes standard polynomials, the results in this work extend the 
recently developed theory on polynomial reproduction to the non-stationary 
context. A significant application of the derived algebraic conditions on the 
subdivision symbols is the construction of new non-stationary subdivision 
schemes with specified reproduction properties. 

Keywords: Non-stationary subdivision, Symbols, Exponential polynomials. Exponential- 
generation, Exponential-reproduction. 
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1 Introduction 

Non-stationary subdivision schemes have proven to be efficient iterative algorithms 
to construct special classes of curves ranging from polynomials and trigonometric 
curves to conic sections or spirals. Aim of this paper is to establish the algebraic 
conditions that fully identify the reproduction properties of a given convergent, 
univariate, binary, non-stationary subdivision scheme. 
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Following the notation in [6], we denote by {a'^'^\ k > 0} the finite sequence of 
real coefficients corresponding to the so called k-level mask of the non- stationary 
subdivision scheme and we define by 

aW(^) = ^afz^ k>0, zeC\{0} (1) 

jez 

the Laurent polynomial whose coefficients are exactly the entries of a^'^^ The poly- 
nomial in ([1]) is commonly known as the k-level symbol of the non-stationary sub- 
division scheme. 

Hereinafter we will denote through {S'jj(fe), A; > 0} the linear subdivision operators 
based on the masks {a^'^^ k > 0}, identifying the refinement process 

f^'^+i) := 5,(.)fW, (5.(.)f := J2^t^2, ff\ k > 0, (2) 

starting from any initial "data" sequence f*^''^ = f := {/j G M, i G Z}. Note 
that, whenever the recursive relation in ([2]) relies on the same mask at each level 
of refinement - namely a*^'^^ = a for all /c > - then the subdivision scheme is said 
stationary. 

Attaching the data f^^^ generated at the k-th step to the parameter values tf^^ 
with 

tfUtSS, and ta-tf) = 2-^ k>0 

we see that the subdivision process generates denser and denser sequences of data so 
that a notion of convergence can be established by taking into account the piecewise 
linear function F^''^ that interpolates the data, namely 

F^'\tl''^) = ft\ GHi, zGZ, k>0, 

where Hi is the space of linear polynomials. If the sequence {F^''\ k >0} converges, 
then we denote its limit by 

gf := lim F^''^ 

and say that gf is the limit function of the subdivision scheme based on the rule ([2]) 
for the data f . Due to the linearity of the refinement rule, we can equivalently check 
convergence of the subdivision scheme when applied to the initial data 6 = {Si^ : 
? G Z} = {. . . , 0, 0, 1, 0, 0, . . . }. If this is convergent with limit = g^, then we have 

gf = - 

jez 

for any initial data sequence f . 

Since most of the properties of a subdivision scheme (e.g. its convergence, its 
smoothness or its support size) do not depend on the choice of the parameter values 

(k) 

tl , these are usually set as 

t\'^:=^, zGZ, k>0. (3) 
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We refer to the choice in Q as to the "standard" parametrization. As it will 
be better clarified later, with respect to the subdivision capability of reproducing 
specific classes of functions, the standard parametrization is not always the optimal 
one. Indeed, the choice 

4'^:=^, leZ, peR, k>0 (4) 

with p suitably set, turns out to be a better selection. In particular, when p = 
we call this parametrization "primal", while in the case p = — | we use the 
term "dual" parametrization, as in For a complete discussion concerning the 
choice of parametrization in the analysis of polynomial reproduction for stationary 
subdivision schemes, we refer the reader to the papers [3l H] . In passing, we recall the 
following definition since it lays the foundations of our results in the non-stationary 
situation. 

Definition 1. A convergent stationary subdivision scheme Sa is generating polyno- 
mials up to degree do if for any polynomial ir of degree d < do there exists some 
initial data q^°^ such that S'^c^^^ = ir. Moreover, q*-*^^ is sampled from a polynomial 
of the same degree and with the same leading coefficient. Additionally, a convergent 
subdivision scheme is reproducing polynomials up to degree dpi, if for any poly- 
nomial TV of degree d < dji and for the initial data tz^^^ = {ir^i), i E Z} it results 

We also recall that a convergent subdivision scheme that reproduces polynomials 
of degree dji has approximation order dR + 1 (see [8] for the proof of this result). 

In this context, the first purpose of this paper is to deal with the extension of the con- 
cepts of polynomial generation and polynomial reproduction to the non-stationary 
setting, leading to the concepts of exponential polynomial generation and expo- 
nential polynomial reproduction. In a few words, a subdivision scheme is said to 
possess the property of reproducing exponential polynomials if, for any initial data 
uniformly sampled from some exponential polynomial function, the scheme yields 
the same function in the limit. The second purpose of this paper is to show that, 
as in the stationary case discussed in [3], the choice of the correct parametrization 
is crucial in the non-stationary setting as well. 

Our analysis brings to algebraic conditions, involving a parametrization of general 
type (jlj), that the symbols of a convergent, binary, non-stationary subdivision scheme 
should fulfill in order to reproduce spaces of exponential polynomials. We remark 
that, as the symbol of the scheme changes from level to level and the parametriza- 
tion plays a crucial role in this kind of study, the proofs of the non-stationary case 
are often more difficult than in the stationary case. This is why much of the results 
previously obtained can not be straightforwardly generalized but require a reformu- 
lation (see also [7]). 

To illustrate the application of the simple but very general algebraic conditions 
derived, we construct novel subdivision symbols defining new non-stationary subdi- 
vision schemes with specified reproduction properties. 
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The outhne of the paper is the following. In Section [2] we start by defining the 
space of exponential polynomials and we characterize some of its special instances 
that turn out to be really useful in applications. Then, in Section [3], we provide 
new theoretical results concerning the algebraic conditions that the symbols of a 
convergent, binary, linear, non- stationary subdivision scheme should satisfy in order 
to reproduce functions from a given space of exponential polynomials. Finally, in 
Section |H novel subdivision schemes are introduced to show how the illustrated 
conditions can be easily applied to work out non-stationary refinement rules with 
desired reproduction properties. 

2 Space of Exponential Polynomials 

Definition 2 (Space of Exponential Polynomials). Let T G Z+ and 7 = (70, 71, ■ ■ ■ , 
7t) with '-jT ^ ^ a finite set of real or imaginary numbers and let the n-th order 
differentiation operator. The space of exponential polynomials Vt,7 is the subspace 

T 

Vt,^ := {/ : M ^ C, / G C^(R) : ^ y,D^ f = 0}. (5) 

0=0 

A characterization of the space Vr,7 is provided by the following: 

Lemma 1. JI]/ Let 7(2;) = Yli'j=Qlj'^^ '^^^ denote by {6'^, r^j^^i ...^^r the set of zeros 
with multiplicity of •y^z) satisfying 

y'-^\ee) = 0, r = 0,--- ,r,-l, ^ = l,---,iV. 

It results 

N 

T = J2re, VT,j:=Span{x'-e''^, r = 0,--- £=!,■■■, N}. 

e=i 

It is worth noting that given the space Vt^j of exponential polynomials, we can 
define the space of exponential splines, say SVr,^, as the space of piecewise functions 
whose pieces are exponentials from the space Vt,^, opportunely joined with a certain 
degree of regularity. 

2.1 Characterization of special Vy/y spaces 

In this section we illustrate some classical examples of function spaces of the kind 
([5]). In all the considered cases we assume n G N and we denote by t a non-negative 
real or pure imaginary constant, namely 

t G {0, s,is \s > 0}. 
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Example 2.1. Let = 3 and 6i = 0, 62 = t, 62, = —t with multiplicities ri = 
n + 1, T2 = = 1. It follows that Tf}^=i,2,3 identifies the set of zeros with 
multiplicities of the polynomial equation 7(z) = 2;"^^ — t^z""*"^ = 0. Since the 
associated differential equation is 

- t^D^'+'fiz) = 0, 

we have 

The importance of this function space lies in the fact that it coincides with the 
following spaces: 

• polynomials linear combination of {l,x, ■ ■ ■ , x", a;""*"^, x""*"^}, whenever t = 0; 

• mixed polynomials-trigonometric functions linear combination of {1, x, ■ ■ ■ , x"", 
cos(sx), sin(sx)}, whenever t = is, s > 0; 

• mixed polynomials-hyperbolic functions linear combination of {l,x, ■■■ ,x"', 
cosh(sx), sinh(sx)}, whenever t = s > 0. 

Example 2.2. We continue by introducing another interesting example of function 
space Vt,7. Let N = 2n + l and 9£,i = 1, ■ ■ ■ , 2n + l the set of values 9i = 0, = 
jt, 92j+i = —jt, j = 1, ■ ■ ■ ,n, with multiphcities ri = 2, = 1, £ = 2, ■ ■ ■ , 2n + 
1. They identify the set of zeros with multiplicities of the polynomial equation 

2n+2 n 

7(2) = 7jZ* = z'^ (z^ — j^t^). Since the associated differential equation is 

i=0 j=l 

2n+2 

J2 l^D'f{z) = 0, 
i=0 

we have 

1^271+2,7 — l-L, a^, e ,e ,...,e ,e |. 

The importance of this function space lies in the fact that it coincides with the space 
of 

• polynomials linear combination of {1, x, ■ ■ ■ , x"^^, x^""*"^}, whenever t = 0; 

• mixed linear polynomials-trigonometric functions linear combination of 
{1, X, cos(sx), sin(sx), ■ ■ ■ , cos(r;,sx), sin(?7,sx)}, whenever t = is, s > 0; 

• mixed linear polynomials-hyperbolic functions linear combination of 

{1, X, cosh(sx), sinh(sx), ■ ■ ■ , cosh(nsx), sinh(nsx)}, whenever t = s > 0. 

Example 2.3. We conclude with the illustration of a last important function space 
of type ([5]). Let = 3, ^1 = 0, 62 = t, 63 = —t with multiplicities Ti = 2, T2 = T3 = 
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n. As T£}£=i,2,3 identifies tlie set of zeros with multiphcities of the polynomial 
equation 7(2) = z^i^z^ — t^)" = 0, associated with the differential equation 

j^i-iy ew'-'^^'f{z) = 0, 

thus 

Vri ™ tx —tx „ tx ^ —tx ^n—l tx „n—l —tx^ 

2„+2,7 = |l,a;, e ,e ,xe ,xe e ,x e |. 

It coincides with the space of 

• polynomials linear combination of {l,x, ■ ■ ■ , x^*^, a;^""*"^}, whenever t = 0; 

• mixed linear polynomials-trigonometric functions linear combination of 

{1, X, cos(sx), sin(sa:), ■ ■ ■ , x""'^ cos(sx), x"""^ sin(sx)}, whenever t = is, s > 0; 

• mixed linear polynomials-hyperbolic functions linear combination of 

{1, X, cosh(sx), sinh(sx), ■ ■ ■ , x""^ cosh(sx), x"""^ sinh(sx)}, whenever t = s > 
0. 

In practice, all the function spaces in the three examples can be used to reproduce 
conic sections. Furthermore, the function space in Example 12.21 has the additional 
advantage of reproducing trigonometric curves of higher order (such as the cardioid, 
the astroid, ...) while that in Example 12.31 has the interesting capability of repro- 
ducing spirals. Obviously, for all of them we can construct the corresponding SVr^'y 
spaces. 



3 Exponential reproducing subdivision schemes 

As it happens for polynomials, functions in the spaces Vt,7 can be generated via 
a subdivision recursion, but at the cost of considering non- stationary subdivision 
schemes. An important example of such non-stationary subdivision schemes is de- 
fined by the following set of symbols 

BL'\z) = 2ll{'-^^^] , k>0. (6) 

i=l V 62'= + ! + 1 / 

As shown in [HI [12], the subdivision scheme with symbols in (E]) generates limit 
functions belonging to the subclass of C^~^ degree-ra L-splines (with n = T — 1 

N 

and T = r^) whose pieces are exponentials of the space Vt,7 where 7(2;) = 

i=l 

YliLi ~ ^iY' ■ These functions, called exponential B-splines, are investigated in 
several papers (see, for instance [9], |Tll |12] and references therein). 

For later use, we note that the symbols in ([6]) for z^''^ := , i = 1, . . . , N, satisfy 
the conditions 

5f(-f) = 0, j;'^ ^ = 0, r = l,...,r,-l. 
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Furthermore, for > 1 it turns out that 

i=iVe2^+i+l/ Ve2fc+i+l/ i=i, i^Ae2^+i + 1/ 

Notice that, when 6i = and ri = + 1, then Bn\z) in ([6]) does not depend on k 
and its symbol is the stationary symbol of a degree-n (polynomial) B-spline 

We now continue with some definitions generalizing [H Definition 1.1, Definition 
4.4.] and successively we extend the conditions in |3] to the binary non-stationary 
situation. Note that what we get is also an extension of [21 Theorem 4.1], since the 
latter one is limited to the context of primal subdivision schemes only. We remind 
the reader that, within this section, the space Vr,^ is the one in Definition |2l 

Definition 3 (E-Generation) . Let {a^'^^z), k>0} be a set of subdivision symbols. 
The subdivision scheme associated with the set of symbols {a^'^^z), k > 0} is said 
to be Vr.'y-generating (for short E-generating) if it is convergent and for f G Vt,7 
and for the initial sequence f(o) := {/(tf ), z G Z}, It results 

lim^,(.)---5,(o)fW = /, feVT,j. 



In the recent paper [TT], the authors proved a crucial result concerning E- 
generation under some mild assumptions on the subdivision symbols. 

Proposition!, fill Theorem 1] Any convergent, non- stationary subdivision scheme 
associated with symbols {a^''\z), k > 0} of finite length that do not have pairs of 
opposite roots, is said to be generating a function in Vp,^ if and only if 



0, r = 0,...,re-l, for ^= e5^, £ = 1, . . . , A^. 



Remark 1. ^45 pointed out in fll^ . Proposition IJ\ states the important fact that a 
convergent, non-stationary subdivision scheme generates exponential polynomials if 
and only if its symbols {a^'^^z), k > 0} are divisible by the symbols of exponential 
B-splines {Bn\z), k > 0} in Moreover, exponential B-splines turn out to be 
the shortest functions satisfying these conditions. 

Next we investigate the E-reproducing capability of a non-stationary subdivision 
scheme. 

Definition 4 (E- Reproduction) . Let {a^'^-'(z), k>Q} be a set of subdivision symbols. 
The subdivision scheme associated with the symbols {a'^^\z)^ k > is said to be 
VT,'7-reproducing (for short E-reproducing) if it is convergent and for f G Vt^^ and 
for the initial sequence 

f{0) := {/(tf)), ieZ}, It results 
lim S'a(fc) ■■■S'a(o)f^°^ = /. 
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Definition 5 (Step-wise E- Reproduction). Let {a^''\z), k > 0} be a set of subdivi- 
sion symbols. The subdivision scheme associated with the symbols {a^^\z), k > 0} 
is said to be step-wise VT,7-reproducing if for f G Vt,7 and for the sequences 
f(fc+i) ieZ},k>0, It results 

^ ^^j,jf(fe) eqmvalently, /(tf +'^) = f^^f^^^ C^) 

jez 

The next Proposition generalizes [H Corollary 4.5] to the non-stationary situation. 
This is actually already considered in [2], though in the context of "primal" subdi- 
vision schemes only. However, since it is in fact independent of the parametrization, 
we just recall the statement and skip the proof. 

Proposition 2. fS^ A convergent, non- stationary subdivision scheme based on the 
symbols {a^'^^z), A; > 0} which is step- wise VT,7-reproducing is also V^.^y-reproducing 
and vice versa. 



We continue with a technical lemma that will be later used to prove the main result 
of this paper. 

Lemma 2. Let {a^''\z), A; > 0} be such that for r = 0, . . . ,T£ — 1 

(.fr(^ + |)^ = 5^ag)(^-J+pr(.f)2^ (8) 



^ / j€Z 

Then, for r = 0, . . . , — 1 



and 



'zi'^Y. 



E4-li(2jr {4^'"-' = ip- ^n4y, E4li(2j + 1)^ (zf^r^^^ = {py 

(11) 

Proof. Expanding (z + and (z — j + pY in ([8]) we arrive at 

.(^)Vr^V = ^ag)(p_,)^(,f)2., , = o,...,r, (12) 



z, 



while expanding (z + | + ^Y (^ ~ J + pY i^i & at 

s 



('"^'T(l + ^)' = E4li(P-Jr(4V". . = 0,,,,,r. (13) 
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Next, after substituting p := | in ( IT2|) . we expand (2p — JY and obtain 

f=0 ^ ^ ^ ^ jgZ ' ' 

Now, since the above equality is possible if and only if 

jez 

we arrive at the first condition in ( ITO!) i.e. , 

{4^y (pY = E 4^(2^T (^f ^)^^ r = 0, . . . , r, - 1. 

Next, expanding (2j + 1)^ and using the first condition in f llOp we obtain the second 
condition in ( ITU]) for r = 0, . . . , — 1: 

jez f=o ^ ^ jez 

To prove the second part of the Lemma, with the substitution p := | we write ( !T3|) 



as 



4"') (p+iy = Il 4?+i(2p - jT (-f V^'+S . = 0, . . . , r. (14) 
Then, if we expand {2p — j)^ we arrive at 



and therefore, since the above equality is possible if and only if 



we get 

2p 



(.f ) " {2p - ly = J2 4li(2jr" (^f )'^'+\ f = 0, . . . , . = 0, . . . , r, - 1 

which coincides with the first condition in fllip : 

(.f - 1)^ = ai%,{2jy {^f'f'^\ r = 0, . . . , r, - 1. 
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Last, expanding (2j + 1)*" and using the just proven first condition in (!TT|) . we obtain 
the second condition in (fTTl) for r = 0, . . . , — 1 



□ 

We proceed with another useful intermediate result. 

Proposition 3. For the polynomial coefficients {pr,si '^r,s}s=o,---,r satisfying 

r r—l r r—l 

^p,,,(2jr :=n(2j-?) J]d.,.(2jr :=n(2j + l-?)' (15) 

s=0 q=0 s=0 q=0 

it holds 

r r—l r r—l 

s=0 i=0 s=0 i=0 

Proof. The proof is trivial since the first equality simply follows by replacing (2j) 
with p in the first of ( IT5l) . while the second one by replacing (2j) with p — 1. □ 

We are now in a position to state the central result of this paper, which is an 
extension of the main theorem in [3] to the binary, non- stationary context. 

Theorem 1. Let z\ := e2'=+i, £ = I,-- - ,A^. A convergent subdivision scheme 
associated with {a^^\z), ^ > 0} reproduces Vt^^ if and only if for each k > it 
holds 



r = 1, . . . ,r£ - 1, 

1=0 



^ "-^ r = 1, . . . ,Te — 1. 



(17) 



Proof. Let us consider the function /(x) = e^*^^, r G {0, . . . , — 1} with ^£ a root 
of multiphcity of jiz). It obviously holds that / G Vt,^. Let {a*^'^)(2), k > 0} 
be an E-reproducing non-stationary subdivision scheme. Then, from ([7]) it holds 
(splitting the even and odd elements) 



i-j+P 
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or, equivalently 



+ 1 + I y = e^5^ <^^S - 3 + 1^)^ {zff'^' . (19) 

Taking ( fTSfTOj) for r = and summing them up, we obtain 

iGZ jez jez 

as well as 

iez j&L jez 



which are the first relations in (fT71) . To prove the other relations in (fT7|) we proceed 
in the following way. For rG{l,...,r^ — 1}, taking the r — th derivative of 



after multiplication by we have 

jeZ \g=0 / j€Z \q=0 

Therefore, using Lemma |2] and Proposition E] for several times, we obtain 



jGZ \g=0 



5^agVi[2j + l-(r-l)]|n(2j + l-?))( 

9=0 



r-1 

s=0 jGZ 



s=o jez 
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r-1 

(fc)N2i+l 



+[1 - ir-i)]j2dr^i,sY.^2U^3n4'^) 

( r-1 r-1 r-l 

vY^Pr-lAPY + iP- l)J2dr-lAP - ly - - 1) 5^Pr.-l,.(p)^ 
s=0 s=0 s=0 

r-1 



+ [l-{r-l)]J2dr-iAp-^) 



=0 

'r— 1 r— 1 



{4' J [p-ir- 1)] $^Pr-l,.(p)^ + Yl dr-lAP - 1) 



,s=0 s=0 



r-1 



Hence, 



i=0 



(20) 



r-1 

p— r 



dz 



2 (.f ni^ 



i=0 



Differently, by replacing with in ( l20l) . we obtain 

y ^^^^^ = (.f fpEPr-i.(p)^ - (p- i)x;4-i.(p- 1)^+ 

V s=0 s=0 
r— 1 r— 1 \ 

-(r-1) Y^pr-iAvY - [1 - (^ - i)]E ^r~iAp - 

s=0 s=0 J 

(r-1 r-1 \ 

Y,Pr-iApy-J2'^r-iAp-'^y] 
s=0 s=0 / 

= 0, 

which concludes the proof of the necessary part of the claim. 

To prove the converse direction similar arguments can be used. □ 

Remark 2. An important consequence of the previous Theorem is that it allows us 
to identify, whenever it exists, the correct parametrization for a given non-stationary 
subdivision scheme. In fact p has to satisfy 

a('^)(zf ) = 2(^f )^ i = l,...,N. (21) 
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~ (k) 

On the other hand, in case for some i it results = 1 and > 1 for all i = 
1, ■ ■ ■ , TV, by fixing r = 1 we get simpler conditions to be satisfied by p: 

p 2 dz ' 2(4")'"' ''^ ' ' ^ ' 

Definition 6 (Symmetric Schemes). A k-level subdivision mask a^^^ is called odd 
symmetric if a^^} = af\ z G Z and even symmetric if a^!^] = af^-^^, z G Z. In terms of 
k-level symbol these conditions read as a'^^^z) = a^''\z~^) and a'^^\z)z = a^''\z~^), 
respectively. 

From Tlieorem [1] tlie following result follows straightforwardly (see also [1]). 

Corollary 1. If a convergent, non-stationary subdivision scheme with symbols {a^''\z) 
k > 0} satisfies conditions for the root zl''^ = 1 (namely reproduces constants) 
and it is odd ( respectively, even ) symmetric, then it reproduces functions from the 
space of exponential polynomials Vr^y with respect to the primal (respectively, dual) 
parametrization. 

We further observe that Theorem [T] applies to interpolatory symbols with the 
parametrization choice p = and thus, in this special case, it gives back the results 
in [2], Theorem 4.1] and in [TJ Theorem 2.3]. 

Corollary 2. An interpolatory scheme reproduces functions from the space of expo- 
nential polynomials Vt^^ with respect to the primal parametrization p = 0. 

Proof. Since interpolatory symbols satisfy the relation a^''\z) + a*^'^)(— z) = 2, by 
taking the first derivative on both sides of this relation it turns out that 

AaW(^)_ AaW(-2) = 0. (23) 
dz dz 

Hence by evaluating f l23|) at an arbitrary zf'^ := e~2fc+i , it follows that ^a'''^^(z^'^'') = 

Taking into account conditions f|T7|) this holds true only in the case 
/= 0. □ 

We close this section by showing that a shift of z'^ in all the symbols {a^^\z), k > 
0} of an E-reproducing subdivision scheme based on the parametrization p, yields 
the replacement p — )■ p + n in conditions f lT7|) . that is an E-reproducing subdivision 
scheme based on the parametrization p + n. See [5] for the analogous result in the 
stationary case. 

Proposition 4. IfSg^(k), k > is a subdivision scheme with symbols {a^'^^z), A; > 0} 
that reproduces a function space Vr,^, then so does the scheme S^^(k), k > with 
symbols b^'^^z) = a^''\z)z'^, k > for any n E Z,. 
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Proof. According to the Leibnitz rule, the r-th derivative of the symbol h^^\z) is 

i=0 ^ ^ j=0 

For z = —zf^ and i = 1, . . . ,N, each addend of the sum in (12^ vanishes because 



_/k) 

= for z = 0, . . . , — 1, i = 1, - ■ ■ ,N hy assumption, and so 

^^^^^^^^ = 0, r = 0,...,n-l, e = l,---,N. 
To complete the proof, we show by induction that 



dz'' 



2(4'))' " 'n(P + ^-^)' r = 0,...,r,-l, (25) 

i=0 

where p is the parametric shift of S^{k). For r = this holds because h^^\z^^'^) = 
a('=)(zf)) (zf = 2 (^''^y^". Now \et < r < Ti - 1. Since when 2 = zf^ the 
first conditions in f lT7|) can be rewritten in a recursive way as 

6/'+iaW(zf ) _ {p - z) a^''\zP) 

dz'+' ~ dz' ' ^-l'---'^' 

exploiting the fact that (^^^) = (^) + (jli)? we find that 



g^^r+i \ i ) dz^ 



-2; 



i=0 ^ ^ j=0 



\ z / dz'^ 

1=0 ^ ^ j=0 



i=0 i=0 



1=0 ^ ^ i=o 



i=0 ^ ^ i=0 
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1=0 ^ ' j=0 



i=0 ^ ^ ^£ j=0 



T.Q^'^^-^i^'') U^'^-Min-r + i) + ip-i)) 



1=0 ^ ^ j=0 

(p + n-r) (rh^^\zf^] 



Hence, if (1251) is true for r, then it is also true for r + 1, so concluding the proof. □ 

4 Application examples 

In this section we consider some examples of new E-generating and E-reproducing 
non- stationary subdivision schemes. The latter ones are obtained by the application 
of Theorem [1] which is also exploited to compute the correct parametrization to be 
associated with these schemes. 

Example 4.1. Primal scheme generating V^,^ = {1, e*^, e~*^}. 
We start by considering the set of symbols {a^^\z)^ ^ > 0} with 

(fc), ^ „3 (^+1)2 ^2^2^(k)^^^ (2+0(^('')+l))z^+2(2(^,W+2)+3v/^^;('-)+l))z+(2+V2(^,W+l)) 

^1 2 2{^,('=)+l) 4(^(fc)+3+2v/5(;;wTT)) 

(26) 

where, for all /c > 0, v^'^'^ is updated through the recurrence formula 



wW = Wi_L with w(-i)>-l. (27) 

To investigate the generation/reproduction properties of this scheme we use Propo- 
sition [T] and Theorem [T] referring to the function space V^,^ = {1, e*^, e~*^'} 
corresponding to the values 6i =0, 62 = t,6^ = —t and ri = 2, r2 = 1, = \ . 
First of all we verify that a^i\z) does not have pairs of opposite roots. Then, with 
some algebra we easily see that for 

zf^ = l, zf^ = e''''^\ zf = e-'l^'^\ tG{0,.,i.|s>0} 

the conditions 



if{-zf)=% £ = 1,2, 3, 
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and the condition 

are satisfied so that the scheme, indeed, generates the space Vt,^. However, such 
a scheme does not reproduce the whole space above but only the subspace {1, x} 
since conditions in Theorem [T] concerning the roots not fulfilled. 

Finally, we observe that such a scheme is continuous since asymptotically equiv- 
alent [5] to the stationary scheme with symbol 

+ A 2 3 1\ {z + l)\ /I 2 3 1 

>^ + 7^+ o) = + 1) o ^ + 7^ + o 



2^ \8 A 8j 23 ' ' 4 

The next example is to provide a primal scheme reproducing the entire space 
Vt,7 = {1, X, e*^, e~^^}. To this purpose, to keep the E-generation of the space 
Vt,^, we multiply the symbol in (126 p by the factor az"'^ + {l—2a)+az and determine 

the value of a such that the required conditions on Z2''\ z^^ are also verified. 

Example 4.2. Primal scheme reproducing Vt,7 = {1, x, e*^, e~*^}. 

Let 



^ w -i^n ^ 2-t;Wv/2(^ 

2 (2;) = a\ '[z) ■ [az + (1 — 2a) + az), with a — 



2^;(fc)(^;(fc) _ l)^2(t;W + !)■ 

(28) 

Clearly, due to the factor o!{\z), the scheme with /c-level symbol 02*^'* (-2) generates 
the function space Vt,7 = {1, a:, e*^, e~*^}. Our intention is now to show that, 

differently from a{\z), the latter one exactly reproduces that space. First of all, by 
exploiting conditions f l22|) we identify that the correct parametrization is given by 
p = and hence the scheme is primal. Further, since conditions 

a?)(-zf) = 0, a^\zf^) = 2, ^=1,2,3, 

and 

-a, i-z, ) = -a, (z, ) = 0, 

are satisfied, the scheme indeed reproduces the space Vt,^ with respect to the primal 
parametrization p = 0. 

We conclude the example by noticing that the non- stationary subdivision scheme 
( 128|) is of class since asymptotically equivalent [5] to the stationary scheme 
with symbol 

23 I 64^ 16^ 32^ 16^ 64^ 

(£±l)i|'^ -L 1 W 5_^5 _ J7 ^4 I 43^3 I 43^2 _ JAr A^l 

22 "T -L; ^^ 128^ 128^ 64^ 64^ 128^ 128/ ' 
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Aim of the forthcoming example is to define a dual scheme that can reproduce 
the same function space Vt^^ = {I7 e*^, e~*^}. 

Example 4.3. Dual scheme reproducing Vr,^ = {1, x, e*^, e~*^}. 

We consider the /c-level symbol 



"3 {Z) - Z 4 2(i,W+l) 



(29) 



depending on the parameter v^'^^ updated at each step through recurrence (l27j) . 
By applying Theorem [1] referring to the function space Vr^'j = {1, x, e^^, e~*^} 
corresponding to the values 61 = 0, 6*2 = t, 6*3 = — t and ri = 2, r2 = 1, = 1, with 
some algebra we easily see that for 

.f) = l, zf = e*/2^"\ zf = e-*/2^"\ te{0,.,i.|.>0} 

the conditions 

ar)(-.f) = 0, af(.f)=2(.f)-|, £=1,2,3, 
as well as the conditions 

are satisfied so that the scheme reproduces the space Vt,^ with respect to the dual 
parametrization p = —\- 

Since the non-stationary subdivision scheme (!29|) is asymptotically equivalent [5] 
to the stationary scheme with symbol 

[z+lf 5 2 9 5\ {z + lf . / 5 4 1 3 13 2 1 



2^ V 8 4 8/ 22 ' ' V 32 4 16 4 32 

it is as well. 

A last example is given, again, to show the potentialities of Theorem [1] in the de- 
sign of novel non-stationary subdivision schemes with prescribed reproduction prop- 
erties. In fact we now introduce a non-stationary approximating scheme able to re- 
produce a more complicated function space like V^,'/ = {1, X, e^^, e~*^, xe^^, xe~^^}. 
Note that a primal (interpolatory) scheme reproducing this space of functions has 
been recently presented in [TU]. On the contrary, we here focus our attention on a 
dual (approximating) scheme reproducing the same space. 
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Example 4.4. Dual scheme reproducing Vt,^ = {1, x, e*^, e xe*^, xe *^}. 
Let us consider the /c-level symbol 

ay{z) = z 



(30) 



with 

/3 = 4 + 5t;('=) + 2) ^/^^ + 8(i;('=))2 + 17t;W + 6, 

^ = _4^llW+i L (2(^;W)3 + io(^;«)2 + 9v(>^) + 2) + (l6(t;W)2 + 23t;('^) + 6) 



+2 (I6(t;(''))'' + 48(i;(''))3 + 70(z;(*^))2 + 4l7;('^') + 6) 
and v^'^^ defined as in f l27p . 

To verify that the non- stationary subdivision scheme with symbols {a^^\z), k > 0} 
reproduces functions from the space Vt,^ = {1, x, e*^, e~^^, xe*^, xe~*^} corre- 
sponding to the values 9i = 0, 02 = t, 9^ = —t and ri = 2, r2 = 2, rs = 2, we need 
to check that, for 

4^) = 1, 4^) = e*/2^-^\ ^f=e"*/2'="\ tG{0,s,i.|s>0}, 

it results 

af(-4^)) = 0, (;.f )=2(^f ^ = 1,2,3, 

and 

Aaf)(-,f)) = o, A«f)(,f)) = 2p(.fr\ £ = 1,2,3. 

By a few computations it is easy to see that this is exactly the case with respect to 
the parametrization p = —\- 

We finally observe that the non- stationary subdivision scheme (1301) is of class 
since asymptotically equivalent [5] to the stationary scheme with symbol 

(^+1)^ (^^i _ 91 3 I 365 2 _ 91^ I ^ 
26 Vl28 32'* 64 32'*"'" 128/ 

I IN ('_63_^7 _ rib_ & _ i73_^5 , J9^^4 , J79^ 3 _ i73_ ^2 _ , 63 \ 

2^ "T -^J U024^ 1024 1024'* 1024 1024 1024 1024 1024/ ' 
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Application examples of the proposed schemes to data uniformly sampled from 
curves of interest in graphical and engineering applications are given in Figures 
[IJIHI Figure [T] shows that all the three E-reproducing schemes can exactly represent 
conic sections if starting from points uniformly sampled with constant spacing a > 
and choosing v^~^^ = |(e*'^ + e~*'^). Differently, Figure |2] focuses on a reproduction 
capability that characterizes the last presented non-stationary scheme. It consists 
in the exact reproduction of 2D spirals, such as the archimedean spiral and the 
circle involute, which are the curves employed in most gear-tooth profiles. Last, in 
Figure [3] we extend the reproduction of spirals to the 3D case, by considering the 
helix and the conical spiral. While the helix can be reconstructed by all the three 
E-reproducing schemes, the conical spiral can be obtained only by the scheme with 
symbol in ( 130|) . 



Figure 1: Reproduction of conic sections by the subdivision algorithms with /c- level 
symbols in (l28l) . ( l29l) and ( l30l) . The dashed line is the polyline connecting the 
starting points uniformly sampled from the conic section with constant spacing. 
The starting parameter v^~^^ is chosen equal to cos(27r/7), ^, 1, cosh(3/5) for the 
circle, ellipse, parabola and hyperbola, respectively. 
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